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This paper presents simulations of the radiative heat flux imparted on the after-body 
of vehicles entering the Martian atmosphere. The radiation is dominated by CO 2 bands 
emitting in the mid-wave infrared spectral region. This mechanism has traditionally not 
been considered in the design of past Mars entry vehicles. However, with recent analysis 
showing that the CO 2 radiation can be greater than convective heating in the wake, and 
with several upcoming and proposed missions to Mars potentially affected, an investigation 
of the impact of this radiation is warranted. The focus of this paper is to provide a 
better understanding of the impact to aerothermal heating predictions and to provide 
comparisons between NASA’s two main radiation codes, NEQAIR and HARA. The tangent 
slab approximation is shown to be overly conservative, by as much as 58%, for most back- 
shell body point locations compared to using a full angular integration method. However, 
due to the complexity of the wake flow, it is also shown that tangent slab does not always 
represent an upper limit for radiative heating. Furthermore, analysis in this paper shows 
that it is not possible to provide a general knock-down factor from the tangent slab results 
to those obtained using the more rigorous full integration method. When the radiative 
heating is accounted for on the after-body, the unmargined total heat flux can be as high 
as 14 W/cm 2 . 


I. Introduction 

T he radiative component of heat flux on the after-body of vehicles entering the Martian atmosphere 
has traditionally been neglected. However, due to theoretical analyses, * 1 simulations, 5 5 5 ’ and experi- 
ments' highlighting that the mid- wave infrared CO 2 radiation can be significant, the impact to flight vehicles 
needs to be quantified. Understanding the magnitude of the radiation and corresponding uncertainties may 
influence future margin policies, 8 and hence the back-shell material selection and thickness. In this work, 
the magnitude of the radiative heat flux has been calculated on the after-body for two geometries from 
successful Mars missions. These same geometries are being flown on the upcoming InSight and Mars 2020 
missions, and their heat shield sizings may be affected by this analysis. InSight is a 2016 Mars lander aimed 
at measuring the seismic activity on the planet’s surface. The structural geometry is identical to Phoenix 
(70 deg sphere/cone, 2.65 m aeroshell diameter). However, InSight will utilize a faster entry than Phoenix 
(but slower than Pathfinder) at 6.3 km/s. The InSight entry is ballistic with a nominal angle of attack (AoA) 
of zero degrees. Mars 2020 will utilize the same geometry as the Mars Science Laboratory (MSL) which suc- 
cessfully landed on the surface of Mars in 2012. Mars 2020 will include the MEDLI2 instrumentation suite 
that will enable, among other things, the measurement of pressure and in-depth heating to the thermal pro- 
tection system. The in-depth temperature measurements can be used to back out the heating encountered 
by the vehicle during entry. Analysis in this work includes a comparison of the flow-fields generated by the 
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Dplr and Laura computational fluid dynamics (CFD) codes and spectral comparisons between NASA’s 
two main predictive radiation tools, Neqair and Hara. A previous analysis 9 showed that tangent slab was 
overly conservative for after-body radiation calculations for Earth entry, and this work will now examine 
cases relevant to Mars entry. 


II. Description of Predictive Codes 

This section will briefly detail the flow- field and radiation codes used in the present analysis. 

II. A. LAURA 

The Langley Aerothermodynamic Upwind Relaxation Algorithm (Laura) is a high fidelity, structured grid 
flow solver, specialized for hypersonic re-entry physics, utilizing state-of-the-art algorithms for CFD sim- 
ulations. ’ Key elements of Laura include Roe’s averaging and Yee’s Symmetric Total Variation 
Diminishing (STVD) formulation of second-order, inviscid flux. Its non-equilibrium real-gas Navier-Stokes 
flow calculations are parallelized. 

II. B. DPLR 

Dplr uses a finite-volume discretization to solve the reacting Navier-Stokes equations for fluids in thermo- 
chemical non-equilibrium on structured grids. It too is parallelized for efficient computing on large clusters. 
While the software was originally designed for steady-state aerothermodynamic analysis of planetary entry 
vehicles, Dplr has evolved over the years to include a broad spectrum of numerical and physical models 
that enable it to accurately simulate most compressible flows. Additional details on Dplr’s capabilities can 
be found in the references 14, 15, 16. 

II.C. HARA 

The High-temperature Aerothermodynamic RAdiation (Hara) model applied in the present study is dis- 
cussed in detail by Johnston et al. 5 A line-by-line approach is used for atoms and optically thick molecules, 
while a smeared band model is used for optically thin molecules. Hara is based on a set of atomic levels and 
lines from the National Institute of Standards and Technology (NIST) 1 and Opacity Project databases. 
The atomic bound- free model is composed of cross sections from the Opacity project’s online TOPbase, 21 
which were curve fit by Johnston. Hara uses tangent slab as the default option for calculating the wall- 
directed radiative heat flux, with an option for running full angular integration for appropriate cases. Unless 
stated otherwise, results presented using Hara in this work are based on flow-fields generated by Laura. 

II. D. NEQAIR 

Non-EQuilibrium AIR (Neqair) is a line-by-line radiation code which computes spontaneous emission, 
absorption and stimulated emission due to transitions between various energy states of chemical species 
along a line-of-sight. Individual electronic transitions are considered for atoms and molecules, with the 
molecular band systems being resolved for each rotational line. Since the report of Whiting et al., 2 ^ numerous 
updates have been incorporated into Neqair, such as: using the latest version of the NIST atomic database 
(version 5.0), 2 using the bound-free cross sections from TOPbase, 21 incorporating the CO 2 database from 
CDSD-4000, 1 parallelization and improvements to the mechanics of QSS. This latest version of Neqair 
is known as vl4.0. As with Hara, Neqair uses the tangent slab approximation as the default option 
for radiation transport. Furthermore, ancillary utilities enable Neqair to perform full angular integration 
calculations for appropriate cases. Unless stated otherwise, results presented using Neqair in this work are 
based on flow- fields generated by Dplr. 

III. Validation of Relevant Models 


III. A. CO2 Radiation 

The Carbon Dioxide Spectroscopic Databank, CDSD-4000, is the most extensive line list for CO 2 presently 
available. CDSD-4000 is intended to be capable of simulating CO 2 spectra at temperatures up to 5000 K. 
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The sheer number of lines in the database make their full inclusion in Neqair or Hara impractical, so 
different approaches have been taken to implement CDSD-4000 within the codes. 

Neqair has incorporated a reduced form of the database by using a pseudo-continuum approach. The 
Neqair CO 2 model contains 876,000 lines and covers all known CO 2 IR band systems, including the bands 
at 2.0, 2.7, and 4.3 /im. A series of high-temperature CO 2 radiation experimental test cases were identified 
to validate the Neqair CO 2 model. 

To maintain Hara’s efficiency, which is particularly useful for the full angular integration computations 
required in the after-body, the CO 2 IR spectrum is tabulated for a range of temperatures and pressures 
externally from the code (the pressure dependence is required for line broadening above 0.1 atm). These 
tables, which are only computed once, are read and interpolated by Hara for actual shock layer radiation 
computations. Because only a single temperature is used when making the tables, it is assumed that the 
two (translational and vibrational) temperatures in the flow-field are nearly equilibrated, which is typically 
the case, even for the after-body. A total of 500,000 spectral points was found to provide sufficient line 
resolution for the range of temperatures and pressures considered, with weakly radiating regions such as 
3350 - 4000 nm ignored. 

Figure 1 shows a comparison of Neqair and Hara calculations, both based on the same Dplr flow-field, 
for a representative back-shell case. The figure shows that the spectral radiance is in very good agreement 
between the two codes, with the integrated values being within 3% of each other. The tangent slab evaluations 
were within 2%. Thus, any differences greater than a few percent between the two codes in the following 
sections are most likely due to differences in either the flow- field used or the angular integration methods. 
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Figure 1. Comparison of CO 2 back-shell radiation between NEQAIR and HARA. 


As CDSD-4000 contains such a large of number of lines and as the conditions relevant to Mars entry cause 
the radiation on the back-shell to be partially black body limited, care must be taken not to under-resolve the 
spectral grid. For example, the 4 fi m band may initially appear to be black body limited (which would mean 
a sparse spectral grid would be acceptable), as shown in Fig. 2(a). However, as shown in higher resolution 
in Fig. 2(b), the spectrum regularly dips down from the black body limit, and capturing these departures 
from the black body is important for evaluating the total heat flux. These features are not apparent when 
the spectrum is smoothed or down-sampled, as in Fig. 1. 
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Figure 2. Spectral radiance highlighting the level of self-absorption for MWIR CO 2 . Self-absorbed results shown in 
green, optically thin results shown in red. 


III.B. Flow- field Comparisons 

The results presented in this section are taken from representative solutions in the wake to highlight the 
differences between the two flow-fields calculated by Dplr and Laura. Figure 3 shows the profiles of the 
translational temperature, vibrational temperature and the number densities of CO and CO2 from Laura 
and Dplr simulations of the same vehicle and trajectory point. Laura predicts a higher temperature for 
most of the shock layer, while Dplr predicts a higher temperature close to the body. Furthermore, Laura 
is predicting a larger concentration of CO2 in the wake. These differences are due to the use of different 
reaction rates for the dissociation of CO2 and CO and the different surface catalycity models implemented 
in the two codes . 25 




Figure 3. DPLR and LAURA results for (a) T, Tv and (b) mole fraction of CO and CO 2 for an InSight back-shell line 
of sight. A distance of 0 on the x-axis represents the body location. 


III.C. Full Angular Integration 

Due to the complexities of the flow found in the wake of a vehicle, and the more complicated geometries on the 
back shell, using a tangent slab approximation is questionable for after-body radiative heating calculations. 
For more accurate estimation of radiative heat flux on the after-body, a full angular integration is required. 
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Full angular integration is achieved by conducting a numerical integration of radiance at a body point with 
respect to solid angle over all possible lines of sight to/from the body point. The method that has been 
implemented with Dplr, to prepare input data for Neqair, begins with the construction of a unit-radius 
spherical surface octant that is triangulated according to a specified number of points distributed uniformly 
along each edge. In general, four such octants, appropriately transformed, define the hemisphere surface 
(Fig. 4) that can be “seen” by the body point. Two octants commonly suffice when the body point is on 
the center line of an axisymmetric body. The hemisphere is centered at the body point and tangential to 
the surface, with the vertices of the combined octant triangulations defining the lines of sight to/from the 
body point. Note that triangulation provides a more uniformly distributed set of lines than does the use of 
discretized lines of latitude and longitude, but no triangulation of an octant surface can be truly uniform, since 
the elements must be right angle triangles at the three octant “corners” but they tend towards equilateral 
triangles elsewhere. Lines of sight within the flow solution volume grid are constructed by intersecting each 
2-point line between body point and triangle vertex with the outer grid boundary. Finally, the required 
inputs to Neqair are evaluated using a multiblock volume grid form of the Alternating Digital Tree (ADT) 
search method, which interpolates flow-field temperatures and species number densities onto the discretized 
hemisphere lines of sight. Neqair is scripted to run in “line of sight” mode (i.e., no tangent slab evaluation 
is made) on all of these lines. The normal component of radiance from all lines is then integrated with 
respect to solid angle to produce the more thorough estimate of radiative heat flux at the body point. 



Figure 4. MSL geometry showing a unit-radius hemisphere triangulation used to generate lines of sight needed for full 
angular integration calculations. 

A similar methodology for conducting full angular integration calculations has also been implemented 
with the Laura code, to prepare input data for Hara. A brief summary of the implementation follows, 
and further details can be found in Mazaheri et al. The Laura full angular integration (referred to as 
“ray tracing” in Ref. 6) consists of a series of simple vector analyses followed by a searching algorithm that 
finds grid cells for which the lines of sight are passing through. It starts with constructing a series of lines 
of sight spanning the 3-D space from selected surface points outward and determining their intersections 
with volume grid cells. The flow- field solution is interpolated onto each of these lines of sight, and the 
corresponding spectral radiance is calculated with Hara. The integrated intensity along each of these lines 
at the selected surface points is then obtained, and the resultant total radiation incident to these surface 
points computed via integration with respect to solid angle. 
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IV. Impact to Flight 


IV. A. Phoenix Geometry 

The Mars Phoenix spacecraft was launched on August 4th of 2007 and landed successfully on May 25th 
of 2008. Phoenix utilized a 70-degree half-angle sphere-cone rigid capsule fore-body geometry. Similar 
shapes were also successfully used to land payloads for the Viking, Pathfinder and MER missions. NASA’s 
upcoming Discovery Program mission, InSight (due to be launched in March 2016), is based on the same 
sphere-cone geometry as Phoenix. 



Figure 5. 2-D cross-section of InSight detailing normal lines of sight evaluated with NEQAIR. 


This section presents fore-body and after-body heating results corresponding to a proposed maximum 
fore-body heating trajectory for the InSight vehicle. Dplr and Neqair have been used to estimate the 
convective and radiative heating for the body point locations shown in Fig. 5. The results of these calcu- 
lations are presented in Fig. 6. Note that for these calculations, Neqair was run in tangent slab mode on 
axisymmetric flow solutions, both of these assumptions will be shown to be questionable for wake flows in 
the next section. This has been the traditional practice now being reconsidered in this work. Figures 6(a) 
and 6(b) show that at the fore-body locations (stagnation point and shoulder), convective heating is almost 
completely dominant. On the after-body, however, several of the locations show the radiative heating to be 
at a similar magnitude to the convective heating (Figs. 6(c) and 6(e)), while in the case of the umbilical 
location, radiative heating is the dominant mechanism (Fig. 6(d)). Even though the heating values are sub- 
stantially smaller on the after-body, the magnitude is important as the TPS on the after-body is designed 
to handle a much smaller load than the fore-body. Additionally, the radiative heat flux peaks later in the 
trajectory than convection. This is explained by the higher density of CO 2 present at the lower temperatures 
encountered later in the trajectory. Figure 7 shows the impact of the tangent slab approximation versus full 
angular integration as predicted by Laura/Hara over the vehicle surface at one point in the trajectory. 
The tangent slab over-prediction is substantial, between 25% - 30% on average and is particularly large at 
sharp corners, such as the cruise ring step or parachute lid. This phenomenon will be discussed in more 
detail in the next section. 
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Figure 6. Convective and radiative heating axisymmetric/tangent slab calculations by DPLR/NEQAIR, for various 
locations on the InSight geometry. 
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Figure 7. Ratio of full angular integration vs tangent slab for Phoenix geometry calculated with LAURA/HARA. 

IV. B. MSL Geometry 

This section presents radiative heating results for the MSL geometry (which will also be used by NASA for 
the upcoming Mars 2020 mission) for two different trajectory points on the peak 3-sigma convective heat 
load trajectory. Figures 8 and 9 show the shock and wake temperature contours for 71.5 s and 87.5 s as 
calculated by Dplr and Laura respectively. The Dplr and Laura flow-fields show qualitative agreement 
for the dominant flow structure. The 71.5 s point corresponds to peak convective heating on the fore-body. 
A tangent slab calculation at two points (labeled 1 and 3) on the back-shell had previously identified 87.5 s 
as the peak radiative flux on the back-shell. Various interactions of expanding shocks can be observed in the 
instantaneous center-line slices of the unsteady 3-D wakes, creating complex flow- fields. The hot and cold 
spots in the wake change in both position and magnitude throughout the trajectory. Given this variability, 
the actual time of peak back-shell heating may differ were the entire trajectory to be evaluated using full 
angular integration. However, this task would be computationally intensive and was not attempted here. 
The figure also shows the locations of the various body points (co-ordinates detailed in Table 1) used in this 
analysis and the locations of the normal lines of sight from these points. This section will provide the results 
of the radiative heat flux calculations for these points as well as detailing the convergence of the calculations 
as functions of both the extent of the CFD grid into the wake and the number of hemisphere lines of sight 
used in full angular integration. 


T, K 



(a) V = 5259 m/s, p = 8.22e-4 kg/m 3 , T = 176 K, a = 17.1° (b) V = 3624 m/s, p = 2.06e-3 kg/m 3 , T = 198 K, a = 17.5° 


Figure 8. DPLR results for flow-field temperature of an MSL pre-flight design trajectory at (a) 71.5 s and (b) 87.5 s. 
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(a) V = 5259 m/s, p = 8.22e-4 kg/m 3 , T = 176 K, a = 17.1° (b) V = 3624 m/s, p = 2.06e-3 kg/m 3 , T = 198 K, a = 17.5° 

Figure 9. LAURA results for flow-field temperature of an MSL pre-flight design trajectory at (a) 71.5 s and (b) 87.5 s. 


Table 1. Co-ordinates for body point locations used in calculations for MSL geometry 


Number 

Name 

X 

y 

z 



Origin: 

Vehicle Nose 

1 

Aft Seal, Leeside 

0.89863 

0 

2.223282 

2 

Mid-Cone 1, Leeside 

1.40096 

0 

1.84612 

3 

Mid-Cone 2, Leeside 

2.11729 

0 

1.05282 

4 

Parachute Closeout Cone 

2.84241 

0 

0 

5 

Mid-Cone 2, Windside 

2.11729 

0 

-1.05282 

6 

Mid-Cone 1, Windside 

1.40096 

0 

-1.84612 

7 

Aft Seal, Windside 

0.89863 

0 

-2.223282 

8 

Stagnation Point 

0.316756 

0 

-1.06864 
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IV.B.l. Length of the Wake 


This section will determine the impact of how far the CFD grid is extended into the wake has on radiative 
heating. Figures 10(a) and 10(b) show that at t = 87.5 s, Dplr and Laura show quite good agreement for 
T, T v and the number densities of CO 2 and CO for the normal line of sight from body point 2. However, 
when the line of sight is compared for body point 3, as shown in Fig. 10(c) and 10(d), even though the 
magnitude is roughly consistent between the two codes, the spatial variation is not matched well due to the 
different reaction rates implemented. Figures 10(c) and 10(d) also show that, even after extending the flow 
20 m in the wake (which corresponds to 5.07 D), the temperature of the flow is still near 2000 K, where CO 2 
will still radiate significantly in the mid-infrared. Figure 11(a) shows how the tangent slab radiative heat 
flux to the vehicle surface is converging as the length of the wake evaluated is increased. The radiance is 
still not completely converged by 20 m, as the radiative heat flux still changed by 3% compared to the value 
at 18 m. Figure 11(b) shows where the radiative heat flux imparted on the back-shell is originating from in 
the wake for body point 3 using tangent slab. The plot shows that the radiation is steadily increasing from 
20 m to 8 m, then begins to plateau, and finally, a small amount of absorption occurs close to the body. 
The linear dependence near the edge of the domain indicates that the radiance is fairly constant despite the 
decreasing CO 2 concentration. The plateau at 8 m indicates that the radiation is becoming optically thick 
and this is responsible for the convergence observed in Fig. 11(a). However, as the tangent slab is based on 
only a single line of sight, this result does not capture the decreasing solid angle at long distances, nor does 
it account for the off-normal flow structure in the wake. 




Distance into the Wake, m 
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(b) CO 2 , CO for Body Point 2 




(c) T, T v for Body Point 3 (d) CO 2 , CO for Body Point 3 


Figure 10. Comparison of DPLR and LAURA flow-fields for t = 87. 5 s. 
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Figure 11. Convergence of tangent slab radiance as a function of distance into the wake of the DPLR solution for Body 
Point 3, t = 87.5 s. 


Comparisons of the tangent slab and full angular integration calculations for two wake lengths are given 
in Table 2. The full angular integration calculations were performed with the higher (1122 line) resolution. 
The results show that most body points have radiative heat fluxes of several watts per square centimeter, 
and can be as high as 6.7 W/cm 2 or 5.2 W/cm 2 when calculated by tangent slab or full angular integration, 
respectively. Only points 3 and 4 show significant variation with the length of the wake simulated. All the 
other body points cross the shock and reach free-stream conditions regardless of the wake length included, 
and so do not depend upon how far the CFD grid is extended into the wake. The normal lines of sight from 
points 3 and 4, however, remain in the wake their entire length so they display significant sensitivity to wake 
length. The full angular integration, however, accounts for the solid angle subtended by the wake boundary 
and is therefore less sensitive to wake length. In this case the degree of convergence is reduced from around 
40% to 10%. Body point 3 was further analyzed with three different grid sizes at the t=71.5 s point, as 
shown in Table 3. These results indicate that a wake distance of approximately 4 diameters is required to 
capture the radiation emitted onto the back-shell using 3-D integration. 

Table 2. t = 87.5 s. Short Wake Grid is 2.23 D and Long Wake Grid is 5.07 D 


Radiative Heating, W/cm 2 


Line of 
Sight 

Tangent Slab (TS) 

Percent 

Short Wake Long Wake Difference 

2.23 D ( s ) 5.07 D ( z ) (TS Z - TS s )/TS z 

Full Angular Integration (FA) 

Percent 

Short Wake Long Wake Difference 

2.23 D ( s ) 5.07 D (*) (FA Z - FA S )/FA* 

Percent 
Difference 
(TS; - FA0/TS, 

1 

6.705 

6.710 

0.075% 

5.016 

5.190 

3.4% 

23% 

2 

4.479 

4.582 

2.3% 

3.360 

3.559 

5.6% 

22% 

3 

2.772 

4.852 

42.9% 

1.803 

2.021 

10.8% 

58% 

4 

0.754 

1.140 

33.9% 

1.034 

1.224 

15.5% 

-7.4% 

5 

1.353 

1.388 

2.5% 

1.208 

1.279 

5.6% 

7.9% 

6 

2.555 

2.540 

-0.59% 

2.134 

2.143 

0.42% 

16% 

7 

5.081 

5.073 

-0.16% 

3.773 

3.782 

0.24% 

25% 

8 

10.018 

10.086 

0.67% 

9.521 

9.562 

0.43% 

5.2% 
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Table 3. t = 71.5 s, BP 3 


Radiative Heating, W/cm 2 

Length 
of Wake 

Tangent slab 
(TS) 

Convergence 

Difference 

Full Angular 
Integration (FA) 

Convergence 

Difference 

Method 

Difference 

2.28 x D 

1.73 

- 

1.18 

- 

31.8% 

4.21 x D 

2.29 

24.5% 

1.43 

17.5% 

37.6% 

4.93 x D 

2.38 

3.8% 

1.45 

1.4% 

39.1% 


IV. B. 2. Tangent Slab vs. Full Angular Integration 

Typically, radiative heat flux is calculated using a tangent slab method due to its ease of implementation, 
quickness and accuracy for locations of dominant heating, such as the stagnation point. However, as this 
section will show, the tangent slab method does not provide good results in calculating back-shell radiation. 
The tangent slab is inaccurate for after-body calculations due to the complexity of the wake flow. This is 
highlighted in Fig. 12 which shows the radiant intensity (i.e., body normal component of radiance) contour 
plotted on the triangulated surface formed by the intersections of the hemisphere lines of sight with the outer 
grid boundary, for the 8 body points. The origin of each plot corresponds to the body-normal line which is 
the basis for the tangent slab calculation. The full angular heat flux is equivalently calculated by integrating 
over the area of the projection plot with respect to solid angle, while tangent slab contours would appear 
circularly symmetric on this plot. It is evident from the non-uniform and non- symmetrical structure of the 
radiance shown in figures 12(a) to 12(j) that the tangent slab approximation is not valid. Figures 12(k) 
and 12(1) show an example of a body point location where tangent slab should be relatively accurate (the 
stagnation point), as the radiance is symmetric and centered on the origin. Comparing each body point pair 
at the two trajectory times of 71.5 and 87.5 s, it is apparent that the angular dependence of the radiance 
is quite different (e.g., Figs. 12(a) to 12(f)) and in some cases the location of the peak radiance has shifted 
(Figs. 12(g) and 1 2 ( j ) ) . This indicates that a constant tangent slab knock-down factor cannot generally be 
applied over an entire trajectory. 

Results from tangent slab calculations using normal lines of sight from each body point, and the full 
angular integration calculations are given in Table 4. The full angular integration calculations were performed 
with two resolution settings for all body points to give an indication of convergence, one using 306 lines to 
define half the hemisphere, and the other using 1122 lines. The CFD grid was extended 5.07 x D (where D 
is the vehicle diameter) into the wake. The results show that most body points show radiative heat flux of 
several watts per square centimeter, and can be as high as 5.83 W/cm 2 or 3.26 W/cm 2 when calculated by 
tangent slab or full angular integration, respectively. Only points 4 (aft symmetry point) and 8 (stagnation 
line) show consistent results between tangent slab and full angular integration. For the stagnation point, the 
tangent slab method over-predicts the full angular integration value by 7.1%, which is consistent with the 
expected over-prediction of within 10%. Point 4, looking directly into the wake from the base of the vehicle, 
presents a nearly symmetrical geometry and as such is well estimated by tangent slab. Though the tangent 
slab approximation is within 10% of the full angular integration, in this case it is actually an under-prediction 
of the heat flux. This is because the wake is colder on the normal line of sight than it is off angle. This 
makes the point that it is not necessarily safe to assume that tangent slab will always be conservative. Body 
points 1, 2, 3, 6 and 7 all show a reduction of approximately 40% for full angular integration compared to 
tangent slab while body point 5 shows a reduction of 23%. These knock downs are significant; for example, 
this corresponds to a reduction of approximately 2.6 W/cm 2 at body point 1. 

Table 4 also shows the results from the two different resolutions used in the full angular integration 
calculations. All body points except for 1 and 6 are within 4% when the results with the two different 
resolutions are compared, suggesting they have been reasonably well converged. For these cases, the lower 
resolution value is always larger than high resolution, even when the tangent slab approximation is not. 
This means that the coarse grid will usually be a conservative estimate. However, it is still possible that 
sufficiently complicated wake structures may produce an exception to this. Without further full angular 
integration calculations at lower, medium or finer resolutions, it is not clear from this table whether body 
points 1 and 6 have converged. As body point 6 shows the worst level of convergence, two extra resolutions 
were run using 90 and 650 lines, as shown in Fig. 13. This figure suggests that a convergence of less than 
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1% is achieved if more than 600 hemisphere lines of sight are used. 

The convergence of the triangulation itself is determined by evaluating the solid angle calculated based on 
the selected points and comparing it to 2i r. For the lower resolution of 306 lines, all body points are generally 
within 5% of predicting half a hemisphere, while the 1122 line grid shows that the body points are within 
1.4%. While these values are suggestive of the accuracy of the integration calculation, the values in Table 4 
show that they are neither an upper nor a lower bound to convergence. The previously stated conclusion 
that a constant tangent slab to full angular integration knock down factor is not possible is also verified 
with the results presented in this table. For 71.5 s, the knock down factor is on average approximately 40%. 
While for the results presented in Table 2, the knock down factor was on average approximately 25%, with 
body point 3 at 58%. A constant factor would be beneficial in the sense that calculations could be run using 
the significantly more time-efficient tangent slab method, then compensated for with the difference from full 
angular integration by a knock-down factor. However, as this is not appropriate, the full angular integration 
method would need to be applied across the entire trajectory to provide an accurate calculation for total 
heat load at any given body point. 

Table 4. Comparison of Tangent Slab and Full Angular Integration methods at t = 71.5 s 


Radiative Heating, W/cm 2 Percent Difference 


Body Point 

Tangent 

Slab 

(TS) 

Full Angular 
Integration: 
coarse 

(FA C ) [Res. 306] 

Full Angular 
Integration: 
fine 

(FA/) [Res. 1122] 

Tangent Slab 
Accuracy 
(TS - FA/)/TS 

Coarse vs Fine Full 
Angular Integration 
(FA f - FA C )/FA f 

1 

5.83 

3.46 

3.26 

44% 

-6.1% 

2 

4.07 

2.51 

2.44 

40% 

-2.9% 

3 

2.2 

1.443 

1.43 

35% 

-0.9% 

4 

0.8 

0.884 

0.879 

-9.9% 

-0.6% 

5 

1.28 

1.002 

0.998 

23% 

-1.4% 

6 

2.24 

1.57 

1.44 

36% 

-9.0% 

7 

3.22 

1.99 

1.92 

40% 

-3.6% 

8 

2.96 

2.78 

2.75 

7.1% 

-1.1% 
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Radiance norm , W/cm 2 -sr 



(a) 71.5 s / BP 1 


(b) 87.5 s / BP 1 




Radiance norm , W/cm 2 -sr 


(c) 71.5 s / BP 2 


(d) 87.5 s / BP 2 



BP 3 


Radiance norm , W/cm 2 -sr 



BP 3 


D ' ,Hi ance norm , W/cm -sr 


(e) 71.5 s / BP 3 


(f) 87.5 s / BP 3 


14 of 19 


American Institute of Aeronautics and Astronautics 



BP 6 


Radiance norm , W/cm 2 -sr 


0.65 

0.60 

0.55 

0.50 

0.45 

0.40 

0.35 

0.30 

0.25 

0.20 

0.15 

0.10 

0.05 

0.00 



Radiance norm , W/cm 2 -sr 


0.65 

0.60 

0.55 

0.50 

0.45 

0.40 

0.35 

0.30 

0.25 

0.20 

0.15 

0.10 

0.05 

0.00 


(g) 71.5 s / BP 6 


(h) 87.5 s / BP 6 
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(k) 71.5 s / BP 8 


(1) 87.5 s / BP 8 


Figure 12. Surface projected radiance for various body point locations at 71.5 s and 87.5 s. 
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♦Radiative Heat Flux, W/cm2 ♦Percent Change 



Figure 13. 


Convergence of radiative heat flux as a function of number of lines of sight for body point 6 at 71.5 s. 


V. Summary of Results and Code Comparison 

This section provides a summary of the heating magnitudes of interest for the InSight trajectory as well 
as the final converged radiative heating for the MSL geometry at both trajectory times analyzed in this 
paper. Table 5 details the peak convective and radiative heat flux as well as the convective and radiative 
components of the total heat load for the Phoenix/InSight geometry. The table summarizes how convective 
heating is clearly dominant on the fore-body, while the radiation component can be approximately equal 
to convective heat flux on the after-body (Shoulder Seal and Aft Centerline) and even dominant (Umbilical 
location). However, it should be noted that these radiative heating values were calculated with tangent 
slab, which has now been shown here to provide an over-prediction compared to the full angular integration 
method. For a more accurate representation of the impact of radiative heating, full angular integration 
calculations are required. This paper has also highlighted that the tangent slab values cannot be scaled by 
a knock-down factor to approximate full angular integration results. 


Table 5. NEQAIR Peak Convective and Radiation for Phoenix geometry 



Body Point Peak Convective, 

W/cm 2 

Peak Radiative 
(Tangent Slab), W/cm 2 

Convective Heat 
Load, J/cm 2 

Radiative Heat Load 
(Tangent Slab), J/cm 2 

Fore-body 

Stag. Pt. 

65.2 

4.36 

2408 

140 


Shoulder 

69.5 

7.50 

2491 

234 

After-body 

Shoulder Seal 

4.2 

4.56 

176 

152 


Umbilical 

0.48 

1.60 

21 

64 


Cruise Ring 

2.45 

1.54 

120 

60 


Parachute 

5.04 

1.55 

220 

58 


Aft Centerline 

5.21 

4.49 

204 

184 


Tables 6 and 7 show the final converged radiative heating values for MSL calculated with Neqair based 
on Dplr flow-fields. The table also shows code-to-code comparisons against results obtained with Hara 
based on Laura flow-fields for both tangent slab and full angular integration. Laminar convective heating 
values calculated by Laura are also provided to give an indication of the significance for radiative heating 
to the total heat flux. The table shows that excellent agreement is obtained between both sets of codes 
when full angular integration is used, with the radiative heating at both times for all body points being 
within 14% (in general, the difference is much less). However, the agreement when the tangent slab method 
is used is significantly worse (difference as high as 33.4% for body point 4 at t = 71.5 s). This is due to the 
differences in the flow-fields along the normal line of sight and also the manner in which the normal line of 
sight is selected. Small differences in this methodology (such as the exact location of the body point used to 
extract the normal line) when applied to a complicated wake structure, as shown in Figs. 12(c) - 12(f), can 
make a large difference in the tangent slab value. Differences along individual lines of sight become averaged 
out once the entire flow-field is considered, however, leading to more consistent full angular integrations. 
Tables 6 and 7 also show that radiative heating is significantly higher than the laminar convective heating 
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value for body points 1 through 6, and therefore highlights its importance in the design of after-body TPS. 


Table 6. Comparing NEQAIR and HARA for MSL Geometry at t = 71.5 s 



Full Angular Integration 


Tangent Slab 


Laminar 

Body 

Radiative Heating, 

Percent 

Radiative Heating, 

Percent 

Convective Heating, 

Point 

w/. 

cm 2 

Difference 

w/ 

cm 2 Difference 

W / cm 2 


Neqair 

Hara 


Neqair 

Hara 


Laura 

1 

3.26 

3.24 

-0.62% 

5.83 

6.25 

- 7 . 2 % 

0.86 

2 

2.44 

2.39 

2.0% 

4.07 

4.06 

0 . 2 % 

0.09 

3 

1.43 

1.52 

-6.3% 

2.2 

2.13 

3 . 2 % 

0.46 

4 

0.88 

0.973 

-10.6% 

0.8 

1.07 

- 33 . 4 % 

0.17 

5 

1.00 

1.06 

-6.0% 

1.28 

1.43 

- 11 . 7 % 

0.40 

6 

1.44 

1.45 

-0.7% 

2.24 

2.22 

0 . 9 % 

1.90 

7 

1.92 

1.83 

4.7% 

3.22 

2.68 

16.8% 

12.15 

8 

2.75 

2.55 

7.3% 

2.96 

2.75 

7.1% 

51.35 



Table 

7. Comparing NEQAIR and HARA for MSL Geometry at t 

= 87.5 s 


Full Angular Integration 


Tangent Slab 


Laminar 

Body 

Radiative Heating, 

Percent 

Radiative Heating, 

Percent 

Convective Heating, 

Point 

W/cm 2 

Difference 

w/ 

cm 2 Difference 

W/cm 2 


Neqair 

Hara 


Neqair 

Hara 


Laura 

1 

5.19 

5.23 

- 0 . 8 % 

6.71 

7.29 

-8.6% 

0.74 

2 

3.56 

3.57 

- 0 . 4 % 

4.58 

4.60 

-0.4% 

0.18 

3 

2.02 

2.14 

- 6 . 0 % 

4.85 

5.73 

-18.0% 

0.26 

4 

1.22 

1.39 

-13.9% 

1.14 

1.46 

-28.4% 

0.08 

5 

1.28 

1.39 

-8.6% 

1.39 

1.60 

-14.9% 

0.46 

6 

2.14 

2.23 

-3.9% 

2.54 

2.72 

-7.0% 

1.97 

7 

3.78 

3.93 

-4.0% 

5.07 

5.4 

-6.2% 

4.38 

8 

9.56 

9.12 

4.6% 

10.09 

10.02 

0.7% 

21.7 


VI. Conclusion 

The paper has highlighted back-shell radiation predictions for planned and actual Mars entry geometries 
in terms of the two main radiation flow- field codes, Dplr/Neqair and Laura/Hara, which are used at 
NASA Ames and Langley, respectively. Detailed heating results are given for both the upcoming Mars 
Insight flight and the recent MSL entry. In the code comparison, it was first shown that the CO 2 radiation 
predictions within Neqair and Hara agreed within 2% when the flow- field input was identical. Next, full 
angular integrations conducted were shown to have excellent agreement, with the radiative heating calculated 
at all body points with both sets of codes agreeing within 14%, and in general, the agreement is much better. 
The convergence of the full angular integration was analyzed with respect to the length that the CFD grid is 
extended into the wake and the number of angular points / lines of sight used in the full angular integration. 
For the MSL case, a wake distance of approximately 4 D and greater than 600 angular points were estimated 
to converge the results within 5%. The variation obtained with tangent slab evaluations showed significantly 
greater difference between the two codes, up to 33.4%. Several of the analyses presented in this paper have 
shown that tangent slab does not provide an accurate indication of after-body heating, and cannot be scaled 
in a general way to approximate the more rigorous full angular integration result. 

Even though the absolute values of the after-body heating are relatively small compared to the fore-body, 
the total heat flux can be significant compared to what the back-shell TPS material is designed to withstand. 
The total heat flux used for design is further increased due to the traditionally large uncertainty on back-shell 
calculations. Therefore, it is clear that understanding after-body heating is important for TPS design. This 
is particularly vexing for upcoming missions like InSight and Mars 2020, as they are built-to-print missions 
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of Phoenix and MSL, respectively, and neither mission considered radiative heating in the original design. 
By including after-body radiative heating with appropriate uncertainty, there is potential that the margined 
heating value would exceed the build-to-print design. The results presented in this paper show that radiative 
heating is significant for both the MSL and Phoenix geometry, and therefore also needs to be considered for 
any future Mars mission. Furthermore, the results show that after-body radiative heating can be larger than 
convective heating, and could lead to an un-margined total heat flux on the after-body up to 14 W/cm 2 . 
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